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Abstract 



The photon distribution function of a discrete series of excitations of squeezed coherent 
states is given explicitly in terms of Hermite polynomials of two variables. The Wigner 
and the coherent-state quasiprobabilities are also presented in closed form through the 
Hermite polynomials and their limiting cases. Expectation values of photon numbers 
and their dispersion are calculated. Some three-dimensional plots of photon distribu- 
tions for different squeezing parameters demonstrating oscillatory behaviour are given. 



PACS number(s): 42.50.Dv, 42.65.Ky 



1. Introduction 



The photon and quadrature statistics of nonclassical states of light such as squeezed 
states 0, ||, even and odd coherent states |3|, || [5], displaced Fock ( or number ) states 
H 0' B or displaced and squeezed Fock states |I0|] , and corresponding multimode 
states as, for example, two- mode squeezed vacuum [11| differs essentially from the 
statistics of light in coherent states [112], [IB], [TJ|, [T5], [RJ which has the Poissonian photon 



distribution and Gaussian quadrature statistics with equal minimal dispersions of the 
both noncorrelated quadratures. The distributions for nonclassical light have frequently 

HI 



|19||. For correlated light [20], this has been found in 



oscillatory character |T 
and this phenomenon takes place for generalized correlated states [22] as well. 

In present work, we concentrate on statistical properties of states studied in [^3], |24 
and the aim of the work is to obtain the explicit analytic expressions in terms of clas- 
sical polynomials and in terms of multivariable Hermite polynomials for the Wigner 
quasiprobability, the coherent-state quasiprobability and the photon distribution func- 
tion. The discussed states are related to general formulations of different types of 
nonclassical states given in p5l |26], E7L |28|. These states form a discrete series of ex- 



citations of squeezed coherent states in an analogous way as the displaced Fock states 
can be considered as discrete excitations of the coherent states. They were introduced 
for the diagonalization of a class of transition operators which make the transition from 
the density operator to the quasiprobabilities of the complete Gaussian class by form- 
ing the trace. For the formulation of orthogonality and completeness relations, it was 
convenient to define these states in a nonnormalized form with a parametrization in 
such a way that in the limit of vanishing squeezing parameter they become identical 
with the displaced Fock states. Other limiting cases of these states are the eigenstates 
of the canonical operators and a discrete series of their excitations which has been not 
considered up to now. The squeezed-state excitations are interesting because they com- 
bine and generalize the properties of well-known important states as squeezed states, 
number states, and coherent states. Such kind of states might appear in the process of 
parametric excitation of the quantum electromagnetic field in resonators with moving 

32] , |33| produces squeezing in 



29, 



walls for which nonstationary Casimir effect 
quadrature components both in coherent states and Fock states. 

We will calculate in the present paper the Wigner and the coherent-state quasiprob- 
ability and furthermore the photon distribution function in terms of the multivariable 
Hermite polynomials f34j . The distributions and quasidistributions were preliminarily 
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given in in a more conventional form. It was shown that the multivariable Hermite 
polynomials are useful to describe the statistics and some other properties of nonclas- 
sical states |35], [!6|, |37|, |39|. One also could realize experimental creation of these 



states using methods of state engineering suggested in []4Cj. 

It should be remarked that there exists nonclassical light with different degree of the 
nonclassicality. For example, the photon statistics of the generalized coherent states 
P5| . and their partial cases like Schrodinger cats || is the same Poissonian statistics 
as for usual coherent states. But the photon quadrature statistics of the generalized 
coherent states is essentially different from the Gaussian statistics of the uncorrelated 
quadratures in the coherent states. Nonlinear Kerr effect and other nonlinear Hamil- 
tonians depending only on photon number operator produce in the process of time 
evolution from the initial coherent states the generalized coherent states and never 
change the photon statistics but influence the quadrature statistics of the light, i.e. 
the quasiprobabilities. From that point of view, the superposition of coherent states 
discussed in [|J are closer to classical then Schrodinger cats of even and odd coherent 
states 0] in which both photon number and quadrature statistics differ from ones in 
the coherent states. As we will point out the quantum statistics of squeezed-states 
excitations is essentially nonclassical in both aspects, i.e. the photon statistics and the 
quadrature statistics are essentially different from these statistics in the classical coher- 
ent states. There may be different understandings of notions of classical light. Usually, 
classical light is considered as coherent-state light and the degree of deviation from the 
classical light may be evaluated as the minimal distance to coherent states [41 1 or, in 



another way, by the minimal parameter s of the s-ordered quasiprobabilities for which 
regions of negativity in the quasidistributions begin to appear |42| . 

The paper is organized as follows. In the next section, we consider the definitions of 
the squeezed-state excitations and calculate the normalization constants of these states 
in closed form. In section 3, we express the general products of the squeezed-state exci- 
tations in terms of multvariable Hermite polynomials. In section 4, we give the explicit 
expressions of the Wigner quasiprobability in terms of the multivariable Hermite poly- 
nomials and discuss some properties of the plots of the quasidistributions. In section 
5, we calculate in closed form the photon distribution of the squeezed-state excitations 
and using plots of the distributions for some sets of parameters we demonstrate the 
oscillatory behaviour of the photon distribution. The mean value of the photon number 
is found explicitly and discussed in section 6 as well as dispersions of quadratures and 
uncertainty product. In two appendices A and B, the properties of the multivariable 
Hermite polynomials are reviewed and some new relations of these polynomials to clas- 
sical polynomials like Legendre, Jacobi, and Gegenbauer polynomials which we used to 
discuss the photon statistics of squeezed-state excitations are derived. 



2 



2. Definition of a set of 
squeezed-state excitations and 
normalization 



A discrete set of squeezed-state excitations |/3, n; 0, ( n = 0, 1, 2, . . . ) has been intro- 
duced in |23], ^4] for the purpose of a diagonal representation of the complete Gaussian 
class of quasiprobabilities. These states possess interesting properties of orthogonality 
and completeness. An appropriate way is to introduce these states in three steps. The 
first step is the introduction of squeezed vacuum states with the complex squeezing pa- 
rameter £ in the nonunitary approach [[43], |3 and in a nonnormalized form according 
to 



10.0:0' = (1 K| 2 )^exp ( - | a t 2 )K)> 



= (i + icr)^EV^^(-cri2m), (2.0 

m =0 v z " L - 

with the following scalar product leading to a nonstandard normalization 



:i + iei 2 xi + icr 
[l - ce 



(0, 0; e|0, 0; C) = ( V J \ " ~ 1 , (0, 0; -C|0, 0; C> = 1- (2.2) 



The next step is the introduction of a discrete set of excitations of the squeezed vacuum 
states according to 

\0,n;Q = -L( a ]~ Ca ) |0,0;C), (n = 0, 1,2, . . . ) • (2.3) 

The last step is the displacement of these states with the unitary displacement operator 
D(/3,j3*) according to 

n; C) = D(f3, (3*) |0, n; C), D(/3, (3*) = exp {(3a) - (3*a). (2.4) 

The states with arbitrary pairs of squeezing parameters ( and — ( and arbitrary dis- 
placement parameters (3 are mutually orthonormalized and obey a completeness relation 
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as follows 

oo 

</3, m; -CI/5, n; C) = <W, £ 1/3, n; C) (/3, n; -C| = /, (2.5) 

n=0 

where J denotes the unity operator of the Fock space. For vanishing squeezing param- 
eter, one obtains from \(3,n; () the displaced Fock states \/3,n) 

|/3, n;0) = D(/3,/F)-=|0) = £>(/?, /T)|n> = |/3,n). (2.6) 
vn! 



The relations (2.5) form the background for the introduction of the states \(3,n;() in 
the described form with a nonstandard normalization. 

The coordinate representation of the states |/3, n; £) has the form ( ^(g; /3, n; () = 
(g|/3,n;C)) 



^(q;P,n; C) 



1 + lc|2 -'W) 



v^VVi-C; n ^(i-0(i + C*)^ v V2 

1+lc|2 ^ex P {-i4if g -^ + ^ 2 + ^-^ ^* 2 ' 



l-C) 2 ^; l-(2hV V2 V "'"V 4 J' 

(2.7) 

where H n (z) denotes the Hermite polynomials. This formula corresponds to the fol- 
lowing connection between the boson operators a and a' and the canonical operators 
Q and P 

a = ^±S, a^^S, ih[a,ai] = [Q,P)=ihI, (2.8) 
v In V In 



and can be calculated using (q\Q = q(q\ and (q\P = —ih-^(q\ and the normalization 
(qW) — $(<1 — q')- The momentum representation of the states \/3, n; () will be given in 
section 3. 

In this section, we recalculate the normalization constant iV n (|C|) of the states 
|/3, n; () given in in form of the series 



N- 2 (\C\) = (0,n;C|0,n;C) 



I + \C\ 2 \ n +l [n/2] n' / \C\ \ 2k 



This series may be either reorganized or recalculated in the form of the integral 



iV- 2 (|C|)= / dq\^q;P,n;()\ 2 , (2.10) 
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which is expressed in terms of Jacobi polynomials P^' k \z) or Legendre polynomials 
P n (z) as follows ( see Appendix A ) 



N- 2 (\C\) 



i+jcj; 
Ni-ici 



p(0,0) 



2 n 



i+JCg 
i-|Cl s 



i + ICI 2 p ^i + ICI 2 



\ l-ICI 2 n Vi-|C| 



f2.11} 



A representation by special Gegenbauer polynomials Cn (z) with upper index 1/2 is also 
possible. Thus, for the normalized states 



|/3, n; C) norm = N n {\(\)\f3,n; (), 
one obtains in coordinate representation 



(2.12) 



ICP 



1 + ICP 



P„ 



i + ICI 2 
n Vi-|CI 2 



i>(q;(3,n; C), 



(2.13) 



where ip(q; /3,n; () can be taken from Eq.(2.7). The normalization constants do not 
depend on the displacement parameter f3 since the unitary displacement operator is 
cancelled when forming the scalar products as can be seen form Eq.(2.4). The normal- 
ization with the normalization constant iV n (|C|) is only possible for \(\ < 1. However, 
the states \f3, n; are well defined also for \(\ > 1 and can be used for auxiliary purposes 
as, for example, for the formulation of completeness relations p3|, E3, El. 
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3. Calculation of the general scalar 
products 



In this section, we calculate the general scalar product of two arbitrary states |a,m;0 
and \/3,n;0. Using Eq.(2.4) in connection with the product of two arbitrary displace- 
ment operators one can reduce this scalar product to the form 



This formula shows that the general scalar product (a, m; n; can be obtained 
from the special scalar product (0,n;£\(3,n;Q by the substitution (3 — > f3 — a and 
multiplication of the result of this substitution by the phase factor exp {l/2(a*(3— a(3*)}. 

According to Eq.(2.7) one can calculate the reduced scalar product by evaluating 
the following integral in coordinate representation 



(a,m;£\P,n; 




(3.1) 



(0,m;£|/3,n; 




(3.2) 



with the following abbreviation for an integral of the Hermite-Gauss type 




(3.3) 



With the substitution 



1 + KI 2 



(3.4) 
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one obtains an integral of the type of Eq.(B.lO) in Appendix B. However, we wrote 
the exponential term in the integrand as a Gaussian with displaced argument with the 
consequence that the result of the integration is a pure polynomial. It can be expressed 
by the two-variable Hermite polynomials as follows 



;i-0(i-g*)far 

i-ce 



H^(yi,y 2 ), 



R 
( 



i-ce 



;i-0(i-e)(i + ici 2 )(i + iei 2 ) 

(i + C*)(i + 



(i-Q(i-e)(i + ici 2 )(i + iei 2 ) 

' \ (i + C)(i+0 

(i-0(e + o 

i + C* 



yi = 



1+0(1 + imp* -cp 
\ 2(1 - o i - ec* ' m ~ 



:i+o(i + ici 2 )/3-^ 

\ 2(i-c) i-ec* 



(3.5) 



This can be written in terms of the usual Hermite polynomials with the following final 
result for the special scalar product in Eq.(3.2) 



(0,m;£|/3,ra; C) 

f (i + iei 2 )(i + iCl 2 ) ' 
I (i-CO 2 ) 

(-i) n / fl+c 



1 r (P + CP*W + CP) 

exp <^ — 



2(1 -co 



V2™+"m!n! \ V 1 - C£* 
{ !^ } (-l)%!n! 
.ft, ./•("' ./)!!" ./)! 



i-ce 



'2,/(l + |eP)(l + |C| 2 )V 



1 + 



\ 2 (i-cO(e + C) 



(/3 + C/?*)Vn- 



i + lCP 



\ 2 (i-ce)(e + o 



(/3* + r/?) J ■ 

(3.6) 



Recall that the general scalar product (a,m;^\/3,n;0 can be obtained from this spe- 
cial scalar product (0,m;£\f3,n;Q by the above mentioned simple substitutions ( see 
Eq.(3.1) ). 

Let us consider special cases of the general scalar product. The special case £ = ( 
is important for the calculation of expectation values for the states \/3, n; () (section 6). 
In this case, one finds from Eq.(3.6) 



(0,m;C|/3,n;C> 
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i + ICI 2 f \£ + (Ff 
^i-|C| 2exp l 2(1 -Id 2 ) 

(-ir ( r~c~ y ~~ 

vm\\Ii-\c\ 2 



i-ICP 



)fi {m,n} 
^ JI(rn -./)!(// -./)!V |u 



-1)- J m!n! 



i + ICP 



m-j 



i + ICI 2 l + Cl* 
\J (i — ICI 2 )C 2 



i + ICI 2 F + C£ 
\J (i - ICI 2 )C* 2 , 



(3.7) 



In case of equal displacement parameter a — /3, the scalar product (/3, m; n; £) 
becomes independent on /3 and can be obtained by setting /3 = in Eq.(3.6), i.e., by 
using the values of the Hermite polynomials for vanishing argument |23J ( H n (0) = 



</3,m;£|An;C> 



^1^1 (-l) fc +'m!n 



' i + ICI s 

i - ce ) 
( £ + c 



Y c+c* 

fe t^ Jfc!Z!*/(m - 2Jfc)!(n - 2/)! V2(l + |£| V V2(l + |C| 2 



5, 



m—2k,n—2l ■ 



(3.8) 

This scalar product is nonvanishing only if \m — n\ is an even number. By setting 
m = n + 2j and m = n + 1 + 2j, one obtains from Eq.(3.7) 



y(i + l£l 2 )(i + ICI 2 ) 



2'\ 



i-ce 

(n + j)! 



^/!(/+j)!(n-2/)!V4(l + |e| 2 )(l + KI 2 ) 

(i + iei 2 )(i + ici 2 )\^ [Re 



2(i + iei 2 ; 

K + CI 2 



e + c vV( n+2 i) !n! 



(n + j)! 



(i-ce) 2 v W i — ce 

e + C y V(n + 2j)!n! / V(l + m(l + |C| 2 ) \ 

2(i + iei 2 )J (n+j)i n v |i-cei /' l ' j 

(/?,n + l + 2j;f|/?,ra;0 = O, (3.10) 

where we have used the Jacobi polynomials P$' k ^ (z) with j = k for the representation 
( see Appendices A and B ). In case of £ = (, one finds from Eqs.(3.9) and (3.10) or 



and 



S 



from (3.7) by setting [3 = 



({ 3,n + 2j; (\{3,n;Q 

i + ICIV c* 



\1-|C| 2 V 1-iCiV (n + j)\ 



n + 2j)\n\ 



pUJ) 



i+Ki! 

i-iCP 



and 



(/3,n + l + 2j;C|An;C) = 0. 



(3.11) 



(3.12) 



The special case j = of this scalar product gives the inverse squared normalization 
constant of the states \{3, n, () which is given in Eqs.(2.9) and (2.11). In case of £ = 
one finds from Eqs.(3.9) and (3.10) the orthogonality relations 



((3,m; -(\(3,n; = 6 m , n , 



(3.13) 



which are the main reason for the choice of the nonstandard normalization of the states 
n;0- 

The states \a, 0; 0) are identical with the coherent states |ft). Therefore, one obtains 
the Bargmann representation of the states \(3, n; () by an analytic function /(«*) as the 
following special case of the general scalar product pi 



f(a*) = exp 



(ft,0;0|/?,n;C> 



l + |C| 2 )iex P |a*/3-|(a*-/3 



' i+jcj: 

2(* 



0* 



(3.14) 



From this relation, one easily finds the coherent-state quasiprobability of the states 
\(3, n; (). We give it only for (3 = because the transition to arbitrary (3 can be made 
from Q(a, ft*) for (3 = by the simple substitutions a — > a — (3 and ft* — > ft* — (3*. 
Thus, one obtains for the states |0, n; () n0 rm taking into account \a, 0; 0) = |a) 



Q(a, ft*) 



1 (a|0,n;C)(0,n;C|a) 



7T 



(0,n;C|0,n;C) 



7T 



exp 



* i C 2 i C * 2 

ftft H ft H — ft 

2 2 



icr 



'i + ICI ; 



ft Fr, 



' i + lCP 

2C* 



ft 



(3.15) 



Figure 1 shows in representation by the canonical variables (q,p) and with H = 1 the 
coherent-state quasiprobability Q(q,p) of the states |0, n; C)norm for the first six values 
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n = 0, 1, . . . , 5 and for the squeezing parameter £ = (3 — V / 5)/2 = 0.381966 from a bird's 
perspective. The plot range is chosen as the maximal range from to 1/ (2tt) which can 
be taken on by the normalized coherent-state quasiprobability Q(q,p). The maximal 
height 1/(2-71") is only reached for coherent states and the square root of the difference 
of the heigth of Q(q,p) to 1/(2-71") for the considered state is a measure for its distance 
to a "classical" state ||1]. Thus, one can see the nonclassicality of states in such kind 
of figures in a very visual way. If one changes the sign of the squeezing parameter (, 
in our case by transition to ( = —0.381966, one obtains a quasiprobability Q(q,p) with 
squeezing axes which are rotated about an angle 7r/2 relative to the primary squeezing 
axes. The seemingly strange value of the squeezing parameter ( = (3— y / 5)/2 was chosen 
for reason that the asymmetry measure (AiV) 3 of the photon number distribution for 
squeezed coherent states takes on the maximal possible negative value that is not very 
relevant for the purposes of this paper. The qualitative features of the pictures for the 
quasidistributions in the central range of the squeezing parameter do not dramatically 
depend on the value of the squeezing parameter. 

The position and momentum representation of the states \/3,n;() can be obtained 
from the general scalar product by using that the states n; () comprise the states \q) 



and \p) as the following special cases |23|, |24| 



9 + ^0;l\ = (2^expg)k), 



2/i'' / v ' r V 2/i 



q + ip 0;-l 



2/j 



: (2/nr)3exp (-*||)b)- 

(3.16) 

In this way, one gets back the position representation of the states \/3, n; () given in 
Eq.(2.7) which was the starting point of our calculations and in an analogous way we 
get the following momentum representation ( ip(p; /3, n; () = (p\/3, n; ) 

V>(p;/3,n;C) 

l + KI 2 / , . [K 




(p+*yf(/3-/n 



v^iVV i+c J n \\ (i+c)(i-c*)fi 



(3.17) 

The displaced Fock states \/3,n) are the special cases of vanishing squeezing parameter 
C of the states \j3, n; 0, i.e. 

|/3, n;0) = D(f3, /3*)-^=|0, 0; 0) = |/?,n), (3.18) 



and, therefore, one obtains the scalar products of displaced Fock states as special cases 
of the general scalar products of the states \/3,n] Q ( see section 5 ). 
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4. Wigner quasiprobability 



It is worth to describe the states \/3,n; () in terms of quasiprobabilities as the Wigner 
quasiprobability |45| or the coherent-state quasiprobability |46], [47] . These quasiproba- 



bilities have been calculated and given without detailed derivations in [2J[ . In particular, 
the Wigner quasiprobability was found there in terms of a finite series including prod- 
ucts of Hermite polynomials. It turns that the series is summed up, i.e. it is possible 
to express the Wigner quasiprobability using the definition 

+oo 

w b>P) = 2f- J ^ ex P^y)^(^-|)^(g + |), (4.1) 

— oo 

by two- variable Hermite polynomials since the integral in Eq.(4.1) is for the considered 
states of the same form as calculated in Appendix B. 

The Wigner quasiprobability for the states \/3, n; () can be obtained from the Wigner 
quasiprobability for the states \0,n;() by a displacement (3 of the distribution in the 
complex plane of the complex variable a that means for the real variables (q, p) by the 
substitutions 



a: 



q + ip 



2h 



a — /3, q 



(4.2) 



Therefore, we restrict us to the calculation of the Wigner quasiprobability for the nor- 
malized states |0,n; Qnorm- From Eq.(4.1) in in view of Eqs.(2.7) and (2.12), one finds 



w( q p) = «( r Ji±cir 4 

{q,P) 2 n n\ \\1-C\J 



i + ICI 2 



cxp 



|(i + C)g + »(i-C)p| s 



2hn 



(4.3) 



with the following abbreviation for an integral of the Hermite-Gauss type 

+oo 



B 



dx H„ 



ICI 2 



n \\ (i-c)(i + on 



X 



9-2) \ H » 



1 + lCI 1 



\ (i + C)(i-C*)» 



X 

" + 2 



exp 



|i-CI 2 ^V2 



(C-C)g + ^|i-Cl 2 P 
1-ICI 2 



(4.4) 
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The exponential function in the integrand is here complemented in a way that it is 
a displaced Gaussian function and the result of the integration is a pure polynomial. 
With the substitution 



i + ICP 



q + 



X 



(4.5) 



^(i + C)(i-C*)^ ' 27' 

one has an integral of the type in Eq.(B.lO) in Appendix B. The result in terms of 
two-variable Hermite polynomials is 



B 



R = 



V2 



-i) n |i -CI 



h 



TV 



v 1-ICI 2 

/ 2C(i-C 



1-ICI 5 



V 



(i + C) , 
li-CKi + ICI 2 ) 

|i-CI 



|i-CI(i + ICI 2 ) \ 
|i + CI 

2C(1-C) 



;i+c*) 



:i + C)(i + ICI 2 )(i-C*)g-^(i + C> 



\ (i-C*)» 



i-ICI 5 



;i + c*)(i + ici 2 ) (i-o?+i(i+c)p 



\ (i-C)ft 



1-ICI 5 



(4.6) 



Thus, expression (4.3) contains no series but only well-known special functions both 
for the normalization constant and for the quadrature-dependent factor. This can be 
expressed also in terms of the series for the usual Hermite polynomials by 

|(l +C )g + i (l_ C ) p |2- 



W(q,p) = — exp 
nir 



(i-icm 

i+ici 2 y 

1-ICI 2 / A (-1)%! / |CI 

A + ICI 2 \ &i' 2 (^-i)'u + ICI 2 



p. 



"vi-ici 5 



1 + ICI 5 



\2C(i-|CI 2 )ft 



((l + C)g + i(l-C)p) 



(4.7) 



or in complex representation 

W(a,a*) : 

f 1 + ICI 2 
I 1-ICI 2 



2 f (a + C«*)(tt* + C*tt) 

^ exp { T=W 



E 



'-lV'ra! 



ICI 



p / l + ICn ^ j! 2 (n-j)!Vl + |C| 5 

u- ICI 2 / 
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i + ICI 2 
^C(i-KI 2 ) 



(a + Ca*) \Hj\ 



1 + IC 



\C*(i-|Cp 



[a 



+ C«) . (4. 



In representation by the canonical variables (q,p) and with h = 1, Fig.2 shows the 
Wigner quasiprobability W(q,p) of the states \0,n;() norm for the first six values n = 
0, 1, ... ,5 and for the squeezing parameter ( = (3 — \/5)/2 = 0.381966 from a frog's 
perspective. The modulus \(\ was chosen in such a way that the relation of the lengths 
of the major to the minor squeezing axes x max and x m i n which are given by 



x j 



\ 



ICI 



i-icr 



x r 



ICI 



\ i+icr 



(4.9) 



becomes equal to a/5 ( another reason for this choice was indicated in section 3 ). As the 
plot range is chosen, the maximal possible range from — 1/7T to 1/tt which can be taken 
on by the normalized Wigner quasiprobability W(q,p). The maximal values +l/ir or 
—1/tt are reached at the coordinate origin q = 0,p = for states with even or odd 
parity, respectively, or in more general cases with displacements for states with definite 
displaced parity correspondingly defined with respect to this displacement. The Wigner 
quasiprobability W(q, p) of the states |0, n; () n orm for the first six values n — 0,1, ... ,5 
and for the squeezing parameter ( = 0.5 from a bird's perspective are given in p3 . 
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5. Photon statistics of the 
squeezed-state excitations 



The photon statistics is described by the photon distribution function p m which is 
expressed in terms of the matrix elements (m\(3,n;() and taking into account \m) = 
|0, m; 0) by 

p w = iV*(|C|)|<0,m;0|/3,n;C>r (5.1) 

The corresponding matrix elements can be taken either in closed form (3.2), (3.5) or as 
series from Eq.(3.6) as the following special case 

(0,m;0|/3,n; 



, ,, |2 a I (WEl (-!) n v? } (-l)'m!n! , r—— 



2c r n -nv 2c 



whereas the normalization constant iV n (|Cl) is given by Eqs.(2.9) and (2.11). 

Let us consider some special cases. For vanishing squeezing parameter (, one obtains 
from \/3, n; C) the displaced Fock states \/3, n; 0) = \/3, n) already in the normalized form 
( see Eq.(2.6) ) and from Eq.(5.2) it follows with the asymptotic expressions for the 
Hermite functions H n (z) — > (2z) n 

< , m;(M = ^(.Bw^ (-y 

V 2 / J m \n\ 7^k V- (m - j )\{n - i ! 
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eXP ("^)f^ ( ^ rXm(l ' 3|2) ' ^(0) = l, (5.3) 



where L v n {z) denotes the associated Laguerre polynomials ( cf., e.g., ). The photon 
statistics of displaced Fock states \(3] n) for n = 0, 1, . . . , 5 and n = 5, 6, . . . , 10 and for 



the squeezing parameter |/3| = y 25/2 = 3.53533 is shown in Figs. 3 and 4. These figures 
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demonstrate the influence of the excitation on the photon distribution structure. As 
we see, for larger excitations the structure becomes more homogeneous. 

In the special case of squeezed coherent states \(3, 0; (), one obtains from Eq.(5.2) 



(0,m;0|/3,0;C) = (1 + ICl') 1 exp { - 



(/3 + CW 



2CY 
V^IV 2 J 



[3 + (P* 



N 2 (\C\) 



i-ICP 



.i + ICI 2 

In particular, one obtains for squeezed vacuum states |0, 0; () 



(5.4) 



(0, 2m; 0|0, 0; C> = (1 + |C| 2 ) 5 ^— r(-C) m , (0, 2m + 1; 0|0, 0; C) = 0, 



2 m m! 



iVo(KI) 



i-ICP 
i + ICP 



(5.5) 



The squeezed vacuum states |O,O;0 are superpositions of even Fock states |2m). The 
states |0, 1; () are squeezed Fock states |1) for which one obtains 



_ 3 J(2m+1)\ 

(0,2m+l;0|0,l;C) = (l + lCl 2 ) l V \ m , (~C) m , (0, 2m; 0|0, 0; C) = 0, 



2 m m! 



n 2 AK\) 



i-ICI 
i + ICI 1 



2 x I 



(5.6) 



These are superpositions of odd Fock states |2m + 1). More generally, one obtains for 
the states |0, n; () 



(0,n + 2m;0|0,n; C) 



n + 2m)!n! W$ 

-cr E 



(n + m)\ 



2 m (n + m)! 



^ /!(/ + m)!(n- 2/)! V4(l + |C| 2 ) 



c + c 



2 m (n + m)! 
(0,n + l + 2m;0|0,n;C) = 0, m>- 

i + ici 2 ^- 1 



> 



n + 1 



iv„ 2 (ICI) 



i-ICI ; 



p 

J. 7I 



i-ICP 



(5.7) 



where [z/] denotes the integer part of v. The states |0, n;Q are only for n = and 
n = 1 squeezed Fock states which can be obtained by applying a unitary squeezing 
operator onto the Fock states |0) and |1), respectively, but they cannot be obtained 
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for n > 2 by applying a unitary squeezing operator to the corresponding Fock states 
\n). This is due to the property that the unitary squeezing operators involves squares 
of the annihilation operator, i.e. a 2 , in the exponent which applied to the Fock states 
|n) with n > 2, do not annihilate these Fock states and the unitary and the nonunitary 
approach become inequivalent. It is, however, possible to determine initial states from 
which |0, n; C) can be obtained by the same unitary squeezing operators as in the cases 
n = and n = 1 from |0) and |1) ( see, || ) or one can first act with equivalent unitary 
or nonunitary squeezing operators onto the vacuum state and have then to take the 
corresponding excitation operators as in the given definition of the states. 

Figures 5, 6, and 7 show the photon statistics of the normalized states \(3,n;Q 
for n = and \/3\ = 5 in dependence on the squeezing parameter ( in case that the 
large axis of the squeezing ellipse is perpendicular to the displacement parameter f3 
( /VvT rea ^ num bers, Figs. 5 and 6 ) and parallel to the displacement parameter f3 ( 
imaginary numbers, Fig. 7 ). The first special case will be called "perpendicular 
squeezing" ( /3/a/C real numbers ) and the second special case "parallel squeezing" 
( fl/VC imaginary numbers ) or we will speak about "perpendicular geometry" and 
"parallel geometry," respectively. Conditionally, one can also speak about "amplitude 
squeezing" and "phase squeezing," respectively, but this is only correct in a certain 
approximation if the quasidistributions are essentially restricted to small sectors of the 
phase plane because amplitude squeezing is squeezing in radial direction and phase 
squeezing is squeezing in angular direction. The characterization of the considered 
special cases by the kind of complex numbers /?/a/C is an invariant characterization 
meaning that it is independent of the chosen position of the coordinate system (q,p) 
in the phase plane [|44 |. In case of perpendicular geometry, we see the oscillations 
of the Fock-state occupation after the main occupation with pulse-like shape for large 
squeezing parameters \(\ ~ 0.75 ... 1. This range of squeezing parameters is amplified in 
Fig. 6. The oscillations can be explained by the neighbourhood to zeroes of the Hermite 
polynomials for real arguments. Figures 5 and 7 show for vanishing squeezing parameter 
( the photon statistics of coherent states ( Poissonian statistics ). Figures 8 and 9 are 
the analogous pictures for n = 10 and displacement parameter = \/l5 = 3.87298. 
These photon statistics also show oscillatory behaviour where in the second case for an 
intermediate range of values of the squeezing parameter \(\ « 0.25. . .0.5 some Fock 
states |n) with low numbers n are occupied and after a gap in the distribution appears 
again a region of occupation with an pulse-like shape. The cooperation of displacement, 
squeezing, and excitation leads to a diversity of possible features in the photon statistics 
that must be studied furthermore. 
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6. Expect ion values for the 
squeezed-state excitations 



The expection value of an arbitrary operator A for the squeezed-state excitations 
\P, n;C) is given by 

(P, n; (\A\P, n; Q (0, n; C[P(A l*)) f ^(A g*)]0, n; Q , , 

(/3,n;OM;C> (0, n; C|0, n; C) ' 1 ' J 

The scalar product in the denominator determines the normalization factor of the states 
n; C) an d is expressed by special Jacobi or Legendre polynomials in Eq.(2.11). A 
possible approach for the calculation of expectation values of ordered powers of the 
annihilation and creation operator, i.e. ordered moments, is the following. First we 
calculate the expectation value of the displacement operator D(a, a*). Due to 

(D(P, P*)) ] D{a, a*)D((3, (3*) = exp(a/3* - a*[3)D(a, a*), (6.2) 

one finds 

D{^a7) = exp(a/r - a*(3) f' n ' X )^ n]< ;l (6.3) 

(0,n;C|0,n;C> 

where the scalar product (0, n; (\a, n; () can be taken from Eq.(3.7) by setting m = n 
and by substituting (3 — > a. For the well-known identities ( e.g., [113] ) 

+s * s f aa* 



exp(aa T ) exp(— a*a) = exp f — — j D(a, a*), 
exp(— a*a) exp(<W) = exp (— D(a, a*), (6.4) 



one has 



° ,v = ( - 1),= {s?w exp (^ )B(a ' a 



a=a*=0 



aa = exp I"— ) ^ ' ' (6 - 5) 



a=ct*=0 
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Hence we find the normally and antinormally ordered moments for the normalized states 

norm 

by 



]i n k 



a 1 'a 



1 

(0,n;C|0,n;C> 

d \ k ( d \ 1 f aa 



da* 
1 



+ exp(^-) (0,n;C|«,n;C) 



a=a*=0 



a a 



(0,n;C|0,n;C) 

d \ h ( d \ l ( aa 



(' 3 * + fc) exp (-^) <o '" ;c|Q '" ;c> L' (6 - 6) 

where we have seperated the displacement part exp(a(3* — a* (3) from the differentiations 
and obtain in such a way by applying the binomial formula the moments as sums over 
powers of the displacement parameters (3 and (3*. However, it is not easy to calculate the 
necessary derivatives in these expressions for arbitrary numbers k and I. They can be 
expressed by multivariable Hermite polynomials but the transition to usual polynomials 
is then very complicated. Therefore, we represent here a second approach which gives 
the possibility to calculate these ordered moments for low order in a more direct way. 
In p3fl , was derived that the operators a(() and (a(— ()) defined by 

a(C) = ^=^, (a( " C))t = ^T^' [0(C) ' (a(_C))t] = Z ' (6 ' 7) 

play the role of the annihilation and creation operator of the states \0,n;() ( see 
Eq.(6.17) in [23| ). With the decomposition 

atO-CK-C)) 1 t K-O^ + CMC) , RR , 

a = / 1 a = , , (0.8) 

one obtains the action of the operators a and a) onto the states |0, n; C) in the form 



ln y-\ VglO, n-l;C)- Cv^TTlO, n + 1; Q 
fl O,n; Q = 

' ' VI + CC* 



tm A v^ + T|Q, n + 1; Q + CV§0, n - 1; Q fao , 

a'\0,n;C) = . , (o.9) 

l > -s/ vi + CC* 



and as a consequence, for example, 
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ata|0,n;C> = ^-{(n - (n + l)(C)\0,n-,0 



-tyj(n + 2)(n + l)\0, n + 2; C) + CVn(n-l)|0, n " 2 5 C>}, 
<m%n;{) = j^— {(n + 1 - n((*)\0, n; C) 

-CV(n + 2)(n + l)|0, n + 2; C) + CVn(n-l)|0, n " 2 - C>}- (6-10) 

With 

(£>(/?, P*))*aD(P, (3*)=a + (31, (D((3, (3*)) ] a ] D(P, /?*) = a f + /3*7, (6.11) 

and with the vanishing of the scalar products of the states |0, n; () and |0, m; () if the 
difference \n — m\ is an odd number, one easily finds 

a = /3, ot = /T. (6.12) 

In the same way, by using Eq.(6.10) and the explicit expressions for the scalar products 
of states |0,n;C) and \0,m;() ( see Eqs.(3.11) and (3.12) ) one finds the expectation 
value of the number operator iV = a) a in the form 

— = (f3,n;C\a^a\(3,n;C} 
{(3,n;(\P,n;() 
= (0,n;(\(al + f3*I)(a + f3I)\0,n;() 
(0,n;C|0,n;C) 

= 1 7TTW\ { ( n - {n + 1)IC|2 ) (™) 

(<" + ^ (iqf ) - (3f) ) } + ^* < 6 - 13 > 

With the following identity for Jacobi polynomials 

(n + 2)P^\z) - nP^{z) = 2{2n + l)P^\z), (6.14) 

which can be proved by using the explicit representations of the Jacobi polynomials 
with equal upper indices ( e.g., Eq.(A.3) in Appendix A ), one obtains from Eq.(6.13) 
the following simple result of the expectation value of the number operator N = a^a 
for the normalized states \(3,n; Qnorm 

y= B + (n + W + | ^ (615) 
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In an analogous way, one obtains from 

o 2 |0,n;C) = Y ^yn(n-l)\0,n-2;C) 

-C(2n + 1)|0, n; () + ( 2 ^/(n + 2)(n + 1)|0, n + 2; ()}, 

a t2 |0,n;C) = yW 2)(n + 1)|0, n + 2; C> 

+C(2n + 1)|0, n; C) + C* 2 - 1) |0, n - 2; C)}, (6.16) 

by using the representation of the scalar products by means of the Jacobi polynomials 
and with the identity in Eq.(6.14) 



- (P,n;(\a 2 \P,n;() 

CL — 



(P,n;C\P,n;C) 
t 

c 



i-ICP 



71-2 



In + 1 + n 



,i-KI 2 J 



v 



p(0,0) 



i+icn 



+ (3\ a t 2 = a 2 *. (6.17) 



The series of functions f n (y) defined by 



fn(y) 



p(0,0) 



(6.18) 



has for the first six values of n the following form 



/o(y) = -(i-y), /i(y) = 0, f 2 (y) 

3(l-y)\l + 3y + y 2 ) 



-,, h{y) 



2(1 - y y 



Hv) 



1 + 16y + 3Qy 2 + 16y 3 + y 4 



, h{y) 



l + 4y + y 2 ' J ° yy ' 1 + % + y 2 ' 
A(l-y) 2 (l + 5y + y 2 ) 



1 + 2Ay + 76y 2 + 24y 3 + y 4 ' 



(6.19) 



and gives a contribution in Eq.(6.17) only for n > 2. The uncertainty of the expectation 
values of the canonical operators Q and P is given by 



% 
2 



|i-CI 2 c+c* 



(AQ) 2 = (Q - Q/) 2 = £ <| (2n + 1)^^ - ^^n/„(|C| 2 ) , 



ft 



(AP) 2 = (P-P/) 2 = - (2n + l) 



i-ICI 2 i-ICP 
i + CI 2 , c + c 



+ 



i-ICI 2 i-IC 



n/ n (IC| 2 ) , (6.20) 
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and the corresponding symmetrized uncertainty correlation by 



AQAP + APAQ = h 



i-ICI 2 



(2n+l + n/ n (|C| 2 )). 



(6.21) 



Now it is easy to calculate from (6.20) the uncertainty sum and the uncertainty product. 
The uncertainty sum ( or total uncertainty ) 



is invariant with respect to rotations of the system, i.e. depends not on the phase of 
squeezing parameter ( but only on its modulus \(\ and increases linearly with the exci- 
tation number n for fixed \C,\. The inequality is simply the relation between arithmetic 
and geometric means. The uncertainty product 



is not rotation-invariant and depends in a complicated manner on the modulus and on 
the phase of the squeezing parameter ( and increases mainly with the excitation number. 
One can analyse Eqs.(6.20) by a computer for low values of the excitation number n 
and finds that each factor (AQ) 2 and (AP) 2 can be reduced below the value h/2 for 
coherent states and reaches for real ( and only real ) positive or negative values of C, in 
the limiting case £ — > ±1, the value zero for (AQ) 2 or (AP) 2 , respectively. Contrary to 
this, the symmetrized uncertainty correlation according to Eq.(6.21) becomes extremal 
for pure imaginary values of the squeezing parameter ( under fixed excitation number 
n. 

In an analogous way, one can calculate the expectation values, for example, of the 
operators a^ 2 a 2 , a 2 a^ 2 or N 2 . However, the expressions obtained are very complicated 
and one cannot find identities for the Jacobi polynomials in such a way that the ex- 
pressions in the numerator and in the denominator shortens in an essential way. We 
do not write down these complicated expressions but it is possible to find the effective 
polynomials in the numerator and denominator for low values of n in the states \/3, n; () 
by using a computer. 




(6.22) 




(6.23) 
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7. Conclusion 



We discussed statistical properties of the squeezed-state excitations both for the quadra- 
ture component of the light and for the photon number. As it turned out the Wigner 
quasiprobability and the photon distribution function for these states are described 
in terms of multivariable Hermite polynomials. The corresponding limiting cases of 
these multivariable Hermite polynomials expressed in terms of Legendre and Laguerre 
polynomials ( and also Gegenbauer and Jacobi polynomials ) gave a possibility both 
to calculate explicitly some series for usual Hermite polynomials and to express the 
normalization constant of the squeezed-state excitations and their scalar product in 
terms of the well-studied classical polynomials. Correspondingly, the coherent-state 
quasiprobability ( Husimi-Kano function ) was expressed in closed form in terms of 
one-variable Hermite polynomials and Legendre polynomials. The influence of dis- 
placement parameter and squeezing parameter onto the photon statistics in the states 
was calculated. 

It turned out that there exists the range of these parameters for which the photon 
distribution function demonstrates highly oscillatory behaviour which is a characteris- 
tics of strongly nonclassical states known for squeezed states, even and odd coherent 
states, and correlated states. The dependence of the average photon number on the 
squeezing and displacement parameter in the squeezed-state excitations was explicitly 
given in terms of Jacobi polynomials. The main results of the work are given in the 
formulae (2.11), (3.15), (4.3), (6.15), and (6.20) in which in closed form are described 
quadrature and photon statistics of the squeezed-state excitations. Using the method 
of constructing the simple superposition states like even and odd coherent states |3j 
( or Schrodinger cat states j| ) it is possible to introduce the superposition states in 
which the partners of the superposition are the excitations of the squeezed coherent 
states introduced in ^4J] and discussed in the present paper. 

One could conclude that, being dependent on the physical parameters like squeezing 
parameter, displacement parameter, and excitation number parameter, the squeezed- 
state excitations demonstrate a reach range of possibilities to influence the photon 
statistics and promise to be realized experimentally in future as other members of the 
nonclassical-state family. 
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Appendix A: 

New polynomials and their relation 
to Jacobi, Gegenbauer, and 
Legendre polynomials 



By comparison of two different methods of calculation of scalar products of the states in 
Eq.(2.12) for the special case of equal displacement parameter we obtained the following 
relation of new polynomials to the Jacobi polynomials P^^\z) ( see, e.g., |34| ) with 



equal upper indices i — j and with transformed argument z as well as multiplied by 
functions of z 

"" 2] (-l)'(n + j)! /xV ,__,„,,,.,/ 1 



ur:^-ww-^y^^h (A - 1) 

where, according to the parity properties of the Jacobi polynomials, both signs of the 
square root \f\ + x are possible, but one has to take the same signs of this square root 
in the argument of the Jacobi polynomials and in the multiplicators. Before proving 
this relation, we make some transformations. With the substitution 



x = - 7T- x 2 , «-> Vl + x = -— , «-> y= — 7 == , A.2 

{i + yr l + y vi + x + i 

one brings Eq.(A.l) into the following form in which we apply it in this paper 

E ,.„ , ( "tf )! M f^)'=(^>f^. (A-3) 



a l\(l+j)l( k n-2l)\\{l + y)V \l + yJ n Kl-y, 
Now, by the substitution 

z - 1 1 + y i a a\ 

z + l 1 — y 

one finds from Eq.(A.3) 

P^(z) - z» T (n + 3)l (—1 (A 5) 

24 



The following relation of the Jacobi polynomials with equal upper indices P^'^(z) 
to the Gegenbauer ( or Ultraspherical ) polynomials C%(z) is known ( |34j], chap. 10.9, 
Eq.(4)) 

(n + 2j)\j\ 

_ (n + j)\(2j)\ (-iffln-fc + j + l) ^^ (A g) 



(n + 2j)\j\V(j + l) k\(n-2k)\ 

where additionally the explicit representation of the Gegenbauer polynomials is inserted. 
There exists another relation between the Jacobi polynomials with equal upper indices 



Pi\z) and the associated Legendre polynomials P£(z) of the following kind [11 



(n + 2j)!VVT 



z* 



However, this relation is uniquely determined only in case of even j = 2m for the two 
possible signs of the square root a/1 — z 2 and one must be cautious when applying this 
relation. The reason is that the associated Legendre polynomials P{ (z) with odd upper 
indices j after the substitution z = cos (9) do not only depend on cos(6 l ) but also depend 
on sin (8). In the particular case j = 0, one obtains from Eqs.(A.6) and (A.7) 

^°' 0) W = Pn{z) = d(z), ( P n {z) = P°(Z) ). (A.8) 

This is one case where the relation to the Legendre polynomials is uniquely defined. 

Let us now prove the given relations. First we prove the transition from the sum rep- 
resentation of Pfj'ti(z) in Eq.(A.5) to the sum representation in Eq.(A.6). By applying 
the binomial formula to (1 — 1/z 2 ) in Eq.(A.5), one finds 

" l; (n + 2j)\^ k< ^ (( + j)!((-fc)!(n-2i)!2 2r y '' 

The interior sum over I can be evaluated as follows 

[n sP (n + 2j)\ _ 2- 2k {2 3 )\T{n-k + 3 + \) 

^ {I + - k)\{n - 2l)\2* j \r(j + l)( n -2k)\ ' l " > 

This can be proved by complete induction from n — >• n + 1 since the result is obviously 
true for n = and all integer j and k > 0. With the decomposition of the factor 
n + 1 + 2j = (n + 1 — 21) + 2(1 + j) in the numerator and the substitution I — > I + 1 
in the second arising sum term, the proof is easily to perform. Thus, it is proved 
that the expression denoted by P^'^(z) in Eq.(A.5) has the representation by the 
Gegenbauer polynomials given in Eq.(A.6) but it is not yet proved by this approach 
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that it is identical with the Jacobi polynomial denoted by the same symbol. However, 
the relations of the Gegenbauer polynomials to special Jacobi polynomials are known 
as already mentioned. Nevertheless, it is instructive to follow a direct transition from 
the explicit general representation of the Jacobi polynomials specialized for equal upper 
indices to the explicit representation by the Gegenbauer polynomials in Eq.(A.6). We 
make this next. 

Starting from the explicit representation of the Jacobi polynomials with equal upper 
indices by applying the binomial formula and the convolution formula for the binomial 
coefficients ( e.g., [f|, Eqs.(6.6) and (6.7) ), one quickly proceeds to the following double 



sum 



a m =o m K n — rn)\(n + j — m)\(j + m)\ 
(n + jY 2 [n £ zn . 2k (2n-2k + 2j)\ 



(n + 2j)!2"^ (n-2k)\ 

{2k,n+j} (—-\Y 

Y~ i L > ca 11) 

r\(2k - r)\(n + j - r)\(n + j - 2k + r)!' V ' ; 



The interior sum can be evaluated in the following way ( see also chap. 4. 2. 5, 
Eq.(29)) 



{ ^ } (-l) r 1 { ^ } (-2) s (2/-s)! 

^ H(2jfc-r)!(Z-r)!(/-2fc + r)! ~ l\(2l - 2k)\ f^ s\(2k - s)\(l - s)\ 



k\l\(l - k)\' 



(A.12) 



Using this sum evaluation with the substitution I = n+j in Eq.(A.ll) together with the 
"doubling" formula for Gamma functions, one obtains P^'^(z) in the form of Eq.(A.6). 
The sum evaluation in Eq.(A.12) can be proved by complete induction from k — > k + 1 
since it is true for k = 0. The sum is transformed in Eq.(A.12) first into another sum. 
This transformation can be made by setting (— l) r = (1 — 2) r and by applying the 
binomial formula to these powers and then after changing the order of summations 
the convolution formula for the binomial coefficients to reduce the double sum to a 
simple sum. It seems that in this transformed form the proof by complete induction is 
easier to perform by multiplying the relation with (2k + 2) (2k + 1) and by using the 
decomposition (2k + 2)(2k + l) = (2k + 2 - s)(2k + 1 - s) + 2s(2A; + 2 - s) + s(s - 1). In 
the general case of the Jacobi polynomials P^'^(z) with unequal upper indices i ^ j, 
one can make analogous expansions in powers of z ( all powers z n ~ l and not only z n ~ 2k 
are present in this case ), where the coefficients are given by finite sums but it was not 
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to see how these sums can be evaluated by a simple expression as it was possible and 
demonstrated for % = j. This is an advantage of the direct derivation. 

Last we consider the relation in Eq.(A.7). The Associated Legendre polynomials 
are defined by f48fl 



-1)', , d l+ i 



Flu , ^L { ^r— { ,-^ 



- | wrBw^w"'^ {AA3) 

With the substitution I = n + j and by applying the doubling formula for the Gamma 
function, one finds 

The comparison of this formula with Eq.(A.6) by applying the doubling formula of the 
Gamma function for (2j)\ reveals the identity in Eq.(A.7) with the discussed problems 
of the sign of the square root. 

Despite the close relations of the polynomials in Eq.(A.l) to Jacobi polynomials 
with equal upper indices, one must look to them as to new original polynomials because 
these relations include nonlinear argument transformations. It seems to be interesting 
to get more relations for these new polynomials as, for example, generating functions, 
recurrence relations, weight factors, and integrals. 
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Appendix B: 

Multivariable Hermite polynomials 



In the Appendix, we review the properties of two-dimensional and multidimensional 
Hermite polynomials following |j35| , |3~6fl , [pCf| . 

The two-dimensional Hermite polynomials H^(y 1 ,y 2 ) are defined by means of the 



generating function [34 



exp{-iafla + ofl»}= £ ^7^ } (?/)- (B.l) 



m,n=0 mW - 

Here a\ and a 2 are arbitrary complex numbers combined into the two-dimensional vector 
a = (ai, a 2 ) : 

2 2 

ai?a = £ diRikdk, aRy = £ aiR ik y k , (B.2) 

i,fc=0 i,fe=0 

and i? is the symmetric 2x2-matrix 

r =(r U r 12 )- ( R3 ) 

V -n-12 -^22 / 



Introducing the notation 

_ -Rl2 



\fRnR22 



(B.4) 



one can obtain |35] the following formula for the two-dimensional Hermite polynomial 
of zero arguments: 

^ } (0, 0) = n mn \{-l)^RTRt {r 2 ~ l) 4 <;:J^ (y^Tf J . (B-5) 
where 

/i mn = min (m,n), (B.6) 

and integers m, n must have the same parity ( otherwise the right-hand side equals zero 
) . Here the Hermite polynomial is expressed in terms of associated Legendre function. 



28 



For coinciding indices, we get 

^ } (0,0) = n!(-deti2)5P B (^==) , (B.7) 

P n (z) being the usual Legendre polynomial. 

Equation (B.5) can be rewritten in the following equivalent forms, using the Jacobi 
polynomials 

or using the Gegenbauer polynomials 



m-j-n 



(B.9) 

For nonzero vector ?/, function H^(yi : y 2 ) can be written as a finite sum of products 
of the usual Hermite polynomials, 



i 



^ / 2R 12 V m!n! / Q \ „ / C2 \ 

(B.10) 

where 

Ci = RuVi + #122/2, C2 = -R122/1 + -R222/2- (B-ll) 

Using Pauli matrix 

**=(j J), (B.12) 

one has for nonzero arguments the relation of the Hermite polynomial to the associated 
Laguerre polynomial 

(n — m+|n — m|) (m — n+|m — n|) 

Hl&'Hyuto) = VmJt—i-lY-yl 5 y 2 a fW. (B.13) 
where 

i/ mn = max(m,n). (B.14) 
For zero argument, we have the formula 

«(o.o) = (ff y t (^g^ (2r) — , (B, 5) 
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where m and n must have the same parity. 

For multivariable Hermite polynomials one can calculate some integrals. Thus, if 
we denote 

m — mi, m2, . . . , rriN, n — ni, n?, . . . tojv, m*, = 0, 1, . . . , (B.16) 
one has 

J dx h£\x)h£ } (Ax + d) eM-xMx + cx) = -j== exp ^-cM^cjH^iy), 

(B.17) 

where the symmetric 2Nx2N-matrix 

R = ( ^ n Ru \ (B 18) 

with NxN-b locks Ru, R22, R12 is expressed in terms of symmetric NxN-matrices 
S, T, M and NxN-matrix A in the form 

R u = S-^SM^S, 
R 22 = T - -TAM^AT, 

# 12 = --TVIM - ^. (B.19) 

Here the matrix A is transposed matrix /l, and R\2 is transposed matrix R 12 . The 
2N-vector y is expressed in terms of N-vectors c and d in the form 

V = iT 1 ( *| ) , (B.20) 

where N-vectors Z\ and 2 2 are 

zi = i(5M- 1 + M- 1 ,S)c 

^2 = ^(TAM^ 1 + M~ x AT)c + Td. (B.21) 
For matrices S — 2, T = 2, the above formula (B.17) yields 

/ dec j J| H mi (xi)H ni ( AikXk + dA\ exp(-xMx + cx) 
•' ^ - : • fc=i 

VdetM 



expfjcM-^Hfflniy), (B.22) 
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with NxN-blocks Ru, R 22 , R12 expressed in terms of NxN-matrices M and A in 
the form 

R u = 2(1 -M- 1 ), 
R 22 = 2(1 -AM' 1 A), 

R 12 = -2AM- 1 . (B.23) 

The 2 N- vector y is expressed in terms of N-vectors c and d in the form of Eq.(B.20) 
with 

z 1 = M- X c, 

z 2 = -(AM" 1 + M~ 1 A)c + 2d, (B.24) 
2 

according to Eq.(B.21). 

In the special case N — 1, the matrices A and M become scalars and one finds from 
Eq.(B.22) 

+00 2 

J dxH m (x)H n (Ax + d)ex V (-Mx 2 + cx) = ^ ex P {jjfi) H ™^ (B.25) 



with 

R = 2 ( 1 " P 



R = z[ 1 / ' *L V (B-26) 



and 

\ / c \ / \ -| (S. l_ 271 d 

2 1\/m \ / f 1 \ 1 / M M 



^2 



ylc 
M 



+ 2^' UJ 2(l-i±f) -2(l-^)d 



(B.27) 



The parameters Ci and £ 2 hi the arguments of the usual Hermite polynomials in the 
representation of the two-variable Hermite polynomials according to Eqs.(B.lO) and 
(B.ll) are then 

Ci = *i, C 2 = ^ 2 . (B.28) 

Thus, we obtain the following explicit representation of this integral by the usual Her- 
mite polynomials 

fM ' ■ 



it GXP V 4M 



^ J dx H m (x)H n (Ax + d) exp(-Mx 2 + cx) 



IM - I \ / /A/ - l 2 \ <g //;!;/! / 2.1 



M V M 



^ j!(m-j')!(n-j')! W(M - 1)(M - A 2 ) 



II- , I , = I //„-; I / C + 2 ^ I • (B 29) 

\2yjM{M-l)) \2^M(M - A 2 ) ) 
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This means that in the considered special case the calculation of zi, z 2 , and R is sufficient 
for the representation of the result of the integral in Eq.(B.25) by sums over the usual 
Hermite polynomials whereas the additional calculation of y\ and y 2 is necessary for 
the representation by the two-variable Hermite polynomials. 
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Figure captions 



Fig.l 

Coherent-state quasiprobability Q{q,p) of the squeezed-state excitations with the pa- 
rameters (3 = 0, C = 0.381966, and n = 0,1,2,3,4,5 corresponding to a),b),c),d),e),f) 
from a bird's perspective ( h — 1 ). 

Fig.2 

Wigner quasiprobability W(q,p) of the squeezed-state excitations with the same pa- 
rameters (3 = 0, C = 0.381966, and n = 0,1,2,3,4,5 corresponding to a),b),c),d),e),f) 
as in Fig.l from a frog's perspective (h — 1). 

Fig.3 

Photon distribution function for (3 = 3.53533, ( = and for excitation numbers 
n — 0,1, ... ,5, i.e., for displaced Fock states. 

Fig.4 

Photon distribution function for f3 = 3.53533, ( = and for excitation numbers 
n — 5, 6, . . . , 10, i.e., for displaced Fock states. 

Fig.5 

Photon distribution function for (3 = 5, n = as the function of the squeezing param- 
eter C, < £ < 1, i.e., for squeezed coherent states and perpendicular geometry ( large 
squeezing axis perpendicular to displacement ). 

Fig. 6 

Photon distribution function for (3 = 5, n = as the function of the squeezing param- 
eter (, 0.75 < ( < 1, i.e., for squeezed coherent states and perpendicular geometry and 
amplification of the range from 0.75 to 1. 

Fig.7 

Photon distribution function for (3 = 5, n = as the function of the squeezing param- 
eter £, < — £ < 1, i.e., for squeezed coherent states and parallel geometry ( large 
squeezing axis parallel to displacement ). 

Fig.8 

Photon distribution function for (3 = 3.87298, n = 10 as the function of the squeezing 
parameter £, < £ < 1, i.e., for squeezed-state excitations with perpendicular geome- 
try. 

Fig.9 

Photon distribution function for f3 = 3.87298, n = 10 as the function of the squeez- 
ing parameter £, < — £ < 1, i.e., for squeezed-state excitations with parallel geometry. 
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